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Abstract 

We study the two-dimensional space of supergravity solutions corresponding to non- 
supersymmetric deformations of the baryonic branch of Klebanov-Strassler. By com- 
bining analytical methods with a numerical survey of the parameter space, we find 
that this solution space includes as limits the softly-broken ]\[ = 1 solutions of Gub- 
ser et al. and those of Dymarsky and Kuperstein. We also identify a one-dimensional 
family of solutions corresponding to a natural non-supersymmetric generalisation of 
Klebanov-Strassler, and one corresponding to the limit in which supersymmetry is 
completely absent, even in the far UV. For almost all of the parameter space we 
find indications that much of the structure of the supersymmetric baryonic branch 
survives. 
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1 Introduction 



One of the primary aims of the study of gauge/gravity duality is to find a dual description 
of realistic field theories such as QCD. This requires generalising the original AdS/CFT cor- 
respondence [I] to cases with less supersymmetry. For example, the Klebanov-Strassler (KS) 
[2] and Chamseddine-Volkov/Maldacena-Nuhez (CVMN) [31 2] backgrounds constitute exact 
globally regular solutions which are dual to J\f = 1 gauge theories. 

The presence of some remaining supersymmetry played a critical role in these successes, 
both in simplifying the search for solutions and in guaranteeing their stability. Despite this, 
considerable progress has been made with respect to the problem of finding dual descriptions 
which completely lack supersymmetry. One natural way that this can be achieved is by 
finding solutions in which a black hole is present, corresponding to a gauge theory at finite 
temperature [5] . See for example (6j [U |9] . 

Alternatively, one can consider field theories in which supersymmetry is softly broken by 
the insertion of relevant operators into the Lagrangian. By using as a starting point theories 
for which the duality is well understood, it is then possible to find dual gravity theories which 
are deformations of the SUSY case, as was achieved, for example, in [9J [TUl EH E2l E3 EH 
[T5l [T6l [T7] . Specifically, the deformed background will match the original one asymptotically 
in the UV. The fact that the deformed backgounds share many of the features, such as 
symmetries, of the SUSY solutions means that the problem of finding solutions is considerably 
simplified. 

In particular this approach was used in [18] to obtain non-SUSY solutions by deforming 
backgrounds on the baryonic branch of KS [19J. Although the SUSY system reduces to a 
single second-order differential equation, in the non-SUSY case it is necessary to solve the 
full equations of motion, consisting of six coupled second-order equations. It was practical 
to solve for asymptotic expansions essentially because of the similarities to the SUSY system 
- the expansions have the same general form. By combining the expansions with numerical 
solutions it was possible to calculate several quantities in the dual field theory, and this 
confirmed that the behaviour was very similar to that in the SUSY case. 

In this paper we establish a more complete understanding of the space of solutions to 
which the solutions of |18j belong. This can be achieved in part as a result of the fact noted 
above, that the non-SUSY solutions share much of the structure of the SUSY baryonic branch. 
By consideration of the asymptotic expansions we find a two-dimensional parameter space 
which includes several previously studied solutions. In addition to the SUSY baryonic branch 
(and its limits, CVMN and KS itself), we also find the non-SUSY solutions of [9] and [T7] as 
limiting cases. By combining the structure described in [S] with that of the SUSY baryonic 
branch, it is possible to describe a generic non-SUSY solution in terms of transitions between 
regions in which SUSY and non-SUSY effects dominate. 

Additionally, we find some interesting special cases, one corresponding to a natural non- 
SUSY generalisation of KS itself, and the other to solutions in which SUSY is no longer softly 
broken and the UV does not match the SUSY case asymptotically. 
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We begin in section [2] by reviewing relevant aspects of the SUSY baryonic branch and 
CVMN solutions. In section [3] we turn to the non-SUSY solutions. First we review the 
solutions of [9], obtained as a deformation of the CVMN background, before moving on to 
the main solutions of interest — the generalisation of the baryonic branch obtained in |18j . 
Section [4] contains the main results of this paper. We first discuss the behaviour of generic 
solutions, and then concentrate on various special cases and limits. Finally we include some 
remarks on aspects of the dual field theory in section [5| 

2 The SUSY system 
2.1 Overview 

Here we present two field theories, which although on the face of it appear different, are in 
fact connected via 'higgsing' (as discussed in [201 [21]). The two theories are firstly that found 
when iV c D5-branes are wrapped on the 2-cycle of the resolved conifold ('theory A'), and 
secondly the baryonic branch of the Klebanov-Strassler quantum field theory ('theory B'). 

Theory A is given by performing a special twisted compactification (to four dimensions) , of 
six dimensional SU(N C ) supersymmetric Yang-Mills with 16 supercharges, preserving only 4 
of them. It was studied in j3[[22j[23] and has a field content (in the four dimensional language) 
consisting of a massless vector multiplet alongside a 'Kaluza-Klein' tower of massive chiral 
and vector multiplets. The form of the Lagrangian, the weakly coupled mass spectrum and 
degeneracies of the theory are written in |22l [23] . The local and global symmetries are 

SU(N C ) x SU{2) L x SU(2) R x U(1) R , (2.1) 

where the R-symmetry is anomalous, breaking U{V) R — > ^2N C - 

Theory B is a quiver with gauge group SU(n + N c ) x SU(n) and bifundamental matter 
multiplets Ai, B a with i, a = 1, 2. The global symmetries are (where again, the R-symmetry 
is anomalous) 

SU(2) L x SU(2) R x U{l) B x U(1) r . (2.2) 

These bifundamentals transform under the local and global symmetries as 

Ai = (n + N c , n, 2, 1, 1, \) , B a = (n + N c , n, 1, 2, -1, \) . (2.3) 

There is a superpotential which can be written as W = ^ejje a/ gtr [AiB a AjBp\. The field 
theory is taken to be close to a strongly coupled fixed point at high energies and it can be 
shown that the anomalous dimension should be ja,b ~ — \- This field theory is known to 
be dual to the Klebanov-Strassler (KS) background [2] and its generalization to the baryonic 
branch [19] . 

The connection between theories A and B is via 'higgsing' as mentioned above. If we 
give a (classical) baryonic vacuum expectation value to the fields (Ai,B a ) and then expand 
around it, we find that the degeneracies and field content of |22l [23] are recovered. 
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In terms of the Type IIB string backgrounds dual to each of the field theories, this weakly 
coupled field theory connection is manifest as a U-duality [20] (and was further studied in 
[2"T] [2U ). The first background (dual to theory A) can be presented using the 

vielbeins 



E x t 


= e 4 dxi, 


E l 


= i e f +9 




2 


E 3 


= -e4 +fe 
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£P = e f + k dp, E e = e^ +h d9, E v = e? +/l sin 0dp, 

(cDi + adfl), E 2 = ]-e^ +9 (uj 2 - asinddp), 



where we have used the following SU(2) left-invariant 1-forms 



(2.4) 



uii = cos t/jdO + sin ifi sin 6d(p, C02 = — sin tpdO + cos ip sin 

ui3 = dip + cos 6 dip . (2-5) 

This means we can write the background and the Ramond-Ramond 3-form compactly as 

i 

F 3 = e "i* [fiE 123 + f 2 E e ^ + h(E 823 + £^ 13 ) + / 4 (^ pie + £ w2 )] , (2.6) 
where we have defined 



E ijk...l =E i AE j AE k A ... A£ l ) 

/! = -2AT c e- fc -^, / 2 = ^e^" 2 V - 2a6 + 1), 

/ 3 = N c e- k ~ h -9(a - b), h = ^e- k - h -%'. (2.7) 

In this setup, the dilaton is a function <&(/}) of the radial coordinate only, and we set a'g s = 1. 
Then the background is written in terms of six functions, {g, h, k, <£, a, b}, which all depend 
on the radial coordinate p only. It is possible to solve the SUSY system using a set of BPS 
equations that can be derived for the above ansatz. 

The family of solutions we will present in section 2.2 correspond to a dual field theory 
deformed by the insertion of an eight-dimensional operator in the Lagrangian which couples 
the field theory to gravity. This calls for a completion in the context of the field theory which 
is achieved on the supergravity side with a U-duality [20J. We will refer to this procedure 
as the 'rotation'. It amounts to a solution generating technique which yields the 'rotated' 
background, in which the vielbeins are 



/-■— e*h idxi, e p = e* +k h*dp, e e = e* +h hid6, = e* +h h± sinOdip, 



1 *_l„ ,„ s r, 1 * 



ei +9 hi(oj 1 + ad9), e 2 = ^-e^ +9 hi (w 2 - asmOdp), 



3 3 = I e f+*^3(& 3 + cos 6dtp). (2i 
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The 'rotation' leaves the RR 3-form invarian iQbut turns on some new fluxes. The new metric, 
RR and NS fields are then 



ds% 



F 3 
F 5 



e 4 



2$ 



/1 3 / 4 

5 

— K- 



[fie 123 + / 2 e^ 3 + h(e e23 + e^ 13 ) + f A {e pie + e plp2 ) 



4> 



fie e<pp - f 2 e pu - f 3 (e 62p + e* lp ) + f 4 (e 183 + e* 23 ) 



ft 3 / 4 

-n^^dt A dx\ A dx2 A dx%, 
h 

2* 



h 



In the above equations we have a new factor defined as 

h=l- K 2 g 2*. 



(2.9) 



(2.10) 



We choose the constant k to be such that the dual QFT will decouple from gravity (corre- 
sponding to careful removal of the eight-dimensional operator) . The choice that allows this is 
k = e~*°°, where <&oo is the asymptotic value of the dilaton for large p. This requirement re- 
stricts us to those solutions in which the dilaton is bounded at large distances. The rationale 
behind this choice is discussed in more detail in [211 EZj - 

2.2 The SUSY solutions 



The background described in (2.4-2.7) results in a system of non-linear, coupled, first-order 



BPS equations (which are derived in the appendix of [2H])- These can be repackaged using a 
certain change of basis functions |29^ 1301 l3l] into a much simpler form where the equations 
decouple: We rewrite the background functions {g,h,k,a,b} in terms of five new functions 
{P, Q, Y, t, a} according to 



4e 



2/1 



P 2 -Q 2 
P cosh t — Q 

sinhT 
P cosh r — Q 



e 29 = Pcoshr-Q, 
N c b = a. 



■2k 



4Y, 



(2.11) 



Then most of the BPS equations can be reduced to algebraic relations between the functions, 
leaving a single decoupled second-order equation for P (referred to in the literature as the 
'master equation'): 



P" + P' 



P' + Q' , P'-Q 



+ 



P-Q P+Q 



4coth(2p- 2p ) 



0, 



(2.12) 



lr The factor of h 3 ^ 4 in (2.9 1 relative to in (2.6 1 simply cancels the factors contained in the new vielbeins 



(2.81. 
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with 



Q = (Q + N c ) coth(2/9 - 2 Po ) + iV c [2pcoth(2/j - 2 Po ) - 1] 



(2.13) 



where Q and p are two integration constants. Each solution to the master equation ( 2.12[ ) 
generically provides us with two backgrounds, related by the U-duality or rotation described 
2lp 



in section 



We will be most interested in the rotated solutions, which correspond to the 
baryonic branch. However, much of what follows will be concerned simply with the behaviour 
of the background functions, and so will apply equally to the unrotated case (corresponding 



to theory A in section 2.1). Additionally, we will at times deal with solutions in which the 



dilaton grows without bound in the UV. Then, as discussed above, we can see from (2.10) 
that we cannot apply the rotation procedure without the warp factor h vanishing]^] 



The master equation (2.12) describes all solutions compatible with the ansatz ( |2. 4 -2.7). 
However, we will restrict our attention to globally regular solutions. In this case we find the 
solutions have an IR (for p — > 0) of the form 

5N r 2N? 



„2/i 



'2k 



hi Ah\ 

— H 

2 15 



hi 



P 2 + 0(p* 



hi 2 



hi 



+ 



4hi 
~ l5~ 
2hi 



l5N r 



h 2 



1 



4iV c 2 
hi 



P 2 + 0( P A ), 



, 16iV c 2 ~, 4n 



9hj 



1 



8N C 
3hi 



P 2 + 0(p 4 ), 



2p 
sinh 2p' 



(2.14) 



where the exact expression for b holds for all p. Aside from the ability to shift the dilaton, 
encoded in 0q, we therefore have a family of solutions parametrised by hi. The second 



integration constant we expect from the second-order equation (2.12) has been fixed to ensure 



regularity. The same requirement also leads us to fix the values of the integration constants 
appearing in (2.12-2.13) as Q Q = —N c and p Q 



0. 



Turning to the UV, we find that for p 

N? /13 



,2 9 



,2/i 



=29 



c + es p + N c (l - 2p) + 



— e 3 p 
4 



c+ 
N 2 



Ap + Ap 2 e-3 p + 0(e 



4(*-*») 



2c4 
3 

1 + 



N 2 

1 _ c 

3c 4 



/25 



20p-8p 2 ) e~l p + 0( e -i p ), 



2p 



sinh 2p 



3N 2 
Ac 2 



2e 



-2p 



2N C 



(1 



3 V 



(2.15) 



+ 0{e-"P), 



2 Or, more generally, a family of backgrounds parametrised by k in ( |2.10[ ). 

3 This does not necessarily mean that we cannot consider these solutions as belonging to the rotated family. 



The issue is in fact slightly more subtle, and we will return to this point in section 2.3 
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with an additional parameter c_ appearing at the next order, giving two non-trivial parame- 



ters. Of course, we require a smooth solution joining the two expansions (2.14-2.15), and this 



can be seen to be the case numerically. However, there is then only one non-trivial indepen- 
dent parameter; given a value for one of {h\,c+,c-}, the requirement that the interpolating 
solution matches both the IR and UV expansions is sufficient to determine the values of the 



other two. This can be seen numerically; a solution found starting from (2.15) with arbitrary 



values of c+ and c_ will generically be singular in the IR, with a divergent Kretschmann 
scalar [2T1. 



2.3 Exploring the baryonic branch 

We saw in section |2.2| that, constrained by the requirement of regularity, and ignoring the 
possible shift of the dilaton, the SUSY solutions form a one dimensional family. It is convenient 



to parametrise the solutions either by hi, which is defined by the IR expansions (2.14), or 



by Of, which is defined by the UV expansions (2.15). The relationship between h± and c+ is 



known only numerically, but for these SUSY solutions we have 



(2.16) 



for large values of c + and hi [24J. As we will see in section [5j in the rotated solutions h\ 
and c+ correspond to the parameter which explores the baryonic branch; we recover the KS 
solution [2] itself in the limit hi,c+ — > oo. We postpone further discussion of this limit until 
section |4~H where we consider its non-SUSY generalisation |17] . 

Taking the opposite limit, c+ — > 0, we find that h\ —> 2N C . This corresponds to the 
Chamseddine-Volkov/Maldacena-Nuhez (CVMN) solution OH]. This is considerably simpler 
than the general case, and exact expressions are known for the functions which describe the 
solution: 



Nr. 



2k 



e 
N 



l, 

sinh 2/}' 



e 2h 

N 

o 4*-4(/>o 



pcoth 2/5 

N 



P* 
sinh 2 2p 



2h sinh 2 2p. 



(2.17) 



Note that while the IR can be obtained simply by setting h\ 



e 

N 



it,. 



:P 4 + 0(p 6 



a 



l--p 2 + 0(p 4 ), e 



4<E>-4</>o 



9' u 11 3 

the UV is qualitatively different from the general case: 

a 2h ^ 



2N C in Q2.14D , 

= l + ^p 2 + 0(/9 4 ), (2.18) 



e 
N 



P 



+ 0(e" 4p ), a = 4pe~ 2p + 0(e 



6p\ 



$ = p + 0(\ogp). 



(2.19) 



Of particular significance here is the fact that the dilaton grows without bound in the UV. 



As anticipated above, this means that we cannot apply the rotation procedure (2.10) 
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In the general case hi > 2N C the system follows the CVMN solution closely, before 
switching to the generic UV fl2l5j ) for large p (figure [l]). That is, we can identify a scale 



Ph x below which (2.17) is almost satisfied, and above which g, h, and k grow exponentially, 
<I> quickly goes to a constant, and a 7^ b. Notice that b is completely unaffected by this; the 
exact result b = 2p/ sinh2/9 holds for all h\. As hi is increased, p^ x moves further into the IR 
(figure[2j). 

We noted above that the rotation procedure could not be applied to the CVMN solution 



(2.17-2.19). Specifically, in section 2.1 we chose a particular value for the constant appearing 
in the warp factor h = 1 — K 2 e 2 *, namely k = e _ *°°. In the CVMN solution <3? grows without 
bound and this identification is no longer possible. Nevertheless, it turns out that there is a 
sense in which we do obtain the (unrotated) CVMN solution by taking the limit hi —> 2N C 
in the (rotated) baryonic branch. To see this, note that as we take the limit hi — > 2N C , we 
find that — > 00, and so k — > 0. In this limit we see that h — )■ 1 and the additional fields 



in (2.9) vanish, returning us to the unrotated system (2.4-2.6) at any finite p. 

More explicitly, in a generic solution on the baryonic branch, the dilaton becomes almost 
constant approximately at the scale ph x (figure [T]) . Provided ph x is large enough that the UV 



expansions are valid, we see from (2.19) that for p < pf ll we have <3? ~ p. Taken together, 
these observations mean that we can write 3>oo ~ Ph x ■ We then find numerically (figure [2]) 
that k 2 ~ e~ 2ph i ~ hi - 2N C -> for hi -> 2N C . 

In effect, taking the limit hi —> 2N C in the rotated solutions simply pushes the scale ph ± to 
infinity, while in the region p < p^ the solution becomes exactly the CVMN one. However, it 
is important to note that the two cases are qualitatively different, and the limit is not entirely 
smooth. In particular, we can expect any quantity which depends on the the UV asymptotics 
of the background to behave discontinuously as we take the limit. 

3 Breaking SUSY 

3.1 Deformation of h t = 2N 2 case 



The CVMN solution [31 H] which we obtain in the limit hi = 2N C (section 2.3) can be 
described in terms of S'0(4) gauged seven-dimensional supergravity. The 50(4) gauge group 
corresponds in the full ten-dimensional description to rotations of the 3-sphere (9,0,1^). 

In order to get a four-dimensional world- volume theory we wrap 5-branes on the S 2 (9, ip). 
There is no covariantly constant spinor on S 2 , so to preserve some supersymmetry we have 
to turn on a gauge field so as to cancel the spin connection of the S 2 in the variation of a 
fermion: 

5* ~ TV = (d, + wJTY* - AjTV)e. (3.1) 

This can be achieved, preserving J\f = 1 SUSY, with an abelian field U(l) C SU(2)l, where 
50(4) ~ SU(2)r x S'C/(2)i|^] In the ten-dimensional description, this corresponds to the 

Alternatively, we could preserve JV = 2 SUSY by choosing the (7(1) to be in a diagonal SU(2)d C 
SU(2) r x 8U(2)l, as in [32[33l[23]. 
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(i) (ii) 




Figure 1: Plots of (i) g (solid) and k (dashed), (ii) e , (iii) log a (solid) and log b (dashed, 
black), and (iv) <£, for the SUSY solutions with 2 < hi < 12, increasing from purple to red. 
Here we set N c = 1 and 4>q = 0. 




' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 logA/ii 
-20 -15 -10 -5 5 

Figure 2: Plot showing the dependence of on A/ii = h\ — 2N C in the SUSY solutions. 
For the purposes of this plot we define ph x by k^p^) = 1/3, corresponding to the transition 
between the CVMN UV (k = constant) and the generic UV (k ~ 2p/3). 



S 



'twist' given by the mixing with the S 2 coordinates 6 and ip in (2.4) 



E 3 ~ ti 3 + cos6dtp. (3.2) 

The resulting solution is singular in the IR. However, we can obtain the regular CVMN solu- 
tion by allowing a non-abelian SU(2) field. In the ten-dimensional description this shows up 



in the additional mixing parametrised by a(p) in (2.4). When a(p) = 1, as occurs for instance 
at the origin in the SUSY solution, the gauge field is pure gauge; the gauge transformation 
which removes the field can be written as a coordinate transformation which removes the 
mixing [331 [33] . 

This solution was generalised in [12] by solving the full equations of motion rather than 
the BPS equations, and by allowing a full 50(4) gauge field. We are interested here in the 
simplest SUSY-breaking deformation of the CVMN solution, where we keep the SU (2) gauge 
group, and introduce a mass term which breaks SUSY. This corresponds to the globally 
regular extremal solutions obtained by Gubser, Tseytlin and Volkov (GTV) [9]. 



For these non-SUSY solutions we no longer have an exact solution as in (2.17), although 
we still have 

2g 2k 

for all p. Instead we must rely on expansions in the IR and UV. In the IR, we have qualitatively 



the same as in (2.18): 



Q = P 2 - (I + V i)p' + 0{P% a = 1 + v 2 p 2 + 0{p% 

e 4$ - 4 ^ = l+(^+^)p 2 + 0(/), (3.4) 
where we have introduced v 2 to parametrise the SUSY-breaking deformation. Comparing to 



(2.18) we see that setting v 2 = —2/3 recovers the SUSY CVMN solution. 

As explained in [9], to obtain a regular UV we need —2 < v 2 < 0. We then obtain 
substantially different behaviour to that in the SUSY case. Adapting the notation of 



— = p + G O0 + 0(-J, a = M a p~ 1 / 2 + 0(p-V 2 ), $ = p + 0(logp), (3.5) 

where the parameters M a and Goo can be considered functions of v 2 - The main qualitative 
difference is the presence of additional terms decaying slower than exponentially in the ex- 
pansions for e 2h and a. This is interpreted in |12j as corresponding to a mass which breaks 
SUSY. 

The effect of the SUSY-breaking deformation is most clearly understood by considering 
a, which is affected at leading order (figure [3] (i)). We see that the non-SUSY solutions are 
characterised by a scale />susy- For p < psiJSY, the qualitative behaviour is that of the SUSY 
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solution, (2.17-2.19), while for p > psusy the non-SUSY UV of (3.5) takes over. For a generic 
non-SUSY solution we can define the deformation to a as 



Aa = a - a S usY- (3.6) 

Then we can think of psusy as the scale at which the deformation Aa, which decays slowly in 
the UV, is of comparable magnitude to asusY, which decays much faster. As a result, psusy 
moves towards the IR as we move further from the SUSY solution (figure [4]). Note that this 
does not relate in a obvious way to the SUSY-breaking scale, which it would be more natural 
to associate with the scale above which Aa has decayed significantly, and which moves into 
the UV as we move further from the SUSY solution. 




Figure 3: (i) Plot of logo against p for —1/10 < Av2 < 1/10, where Av2 = V2 + 2/3, showing 
the transition between a ~ e~ 2p and a ~ ^fp at p ~ /?susy- The dashed curves correspond to 
V2 < —2/3, for which a = at psusy- 

(ii) Plot of a against p for the full range —2 < V2 < 0. Again, the dashed curves correspond 
to V2 < —2/3, for which a has at least one zero. The additional oscillations which are in fact 
present in the case V2 = —2 (purple) are not visible at this scale. 



For V2 > 0, a is always positive, and for V2 = 0, a = 1 for all p. As noted above, this 
means that the gauge field is pure gauge, and we can remove the mixing between the spheres 
by a change of coordinates. Thus in this case the internal geometry is simply S 2 x S 3 . Aside 
from the behavior of a, the UV is otherwise unchanged — the other functions h and still 



behave according to (3.5). 

For V2 < —2/3, a has at least one zero. As V2 is reduced, a picks up more oscillations, 
and in the limiting case there are infinitely many zeros. In this limit the UV of the other 
functions no longer that of (3.5) (see figure [5]). Instead the system approaches the 'special 
Abelian solution' of [9]; 

e 2h 1 

— $^V2p. (3.7) 
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Figure 4: Plots showing the dependence of psusy on Av2 = t>2 + 2/3. For the purposes of 
this plot we define psusy by |Aa| = asusY- The solid blue curve corresponds to V2 > —2/3 
and the dashed red curve to V2 < —2/3. 




Figure 5: Plots of (i) <3? and (ii) e against p, for —2 < V2 < —2/3 (dashed curves) and 
—2/3 < V2 < (solid curves), showing the difference between the generic UV (3.5) and the 
limiting case (|3.7l). 
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3.2 Deformation of general case 



We now turn to the solution which was originally presented in [18J. The aim there was to 
find a non-supersymmetric generalisation, preserving the symmetries and structure of the 
baryonic branch solutions with hi ^ 2N C discussed in section [2} This is analogous to the way 



in which the GTV solutions (section 3.1) generalise the CVMN solution (2.17). 



In the non-SUSY case we can no longer make use of the master equation (2.12), which 



was derived from the BPS equations. Instead we have to solve the full Einstein, Maxwell, 
dilaton and Bianchi equations of the system. This amounts to a system of six coupled non- 
linear second-order equations, together with a first-order Hamiltonian constraint. These are 
included in appendix HI We look for solutions to these equations in the form of IR and UV 



expansions with similar forms to the SUSY case (2.14-2.15). 



In the IR we simply impose that the solution is regular, and that the 2-sphere shrinks to 
zero radius at p = 0, as in (2.14). We then have expansions of the form 



=2s 



4$ 



n=0 

oo 

ra=0 



,,2/i 



n=2 

oo 



„2fc 



n=0 



VnP 



(3.8) 



n=0 



n=0 



Substituting into the equations of motion (A.2-A.7) we find five independent parameters, 
which we take to be ko, fo, &2> V2 and W2- We relabel ko = hi/ 2 and /o = e 4flio , so that we 



can recover the SUSY solution (2.14) by setting 



k 2 



2hi 8N t 



V2 



2 
3' 



8N C 
W2 = — 2. 



5 hi 3 3hi 

After the relabeling, the five independent parameters are^] 



hi 



ko 



V2, 



W2, 



(3.9) 



(3.10) 



and the expansions are qualitatively the same as the SUSY case (2.14 

4Y 2 



J2k 



h\ hi 
— + — 
2 2 

hi 
2 P 



1 



hi 
6 



hi h\ 

2k 2 4A^ C 2 
~~h~i 



Shi 



+ . o + 



hi 



+ k 2 p 2 + 0(p 4 ), 



hi 

,<S>-4>a _ ^ i 



4 



2] P 4 + 0(p 
4 



hi 



+ V Z 2 P 2 + 0( P *), 



l + W 2 p 2 + 0(p 4 ), 



1 + V 2 p 2 + 0(p 4 



(3.11) 



^Notice that hi does not refer to the coefficient of p in the expansion for e , as would be expected from 



the form of (3.8 1. This unfortunate notation should not cause confusion because that term will always be zero 



to ensure regularity. 

6 More complete expressions, both for the IR expansions here and the UV ( 3.16 1, can be found in an appendix 
of EES]. 
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In the UV we use a particular generalisation of the SUSY solutions (2.15): 

oo m oo m 



m=0 n=0 

oo m 



'2k 



m=0 n=0 

oo m 



m)p/3 



m=l n=0 



e2A = EE^ v(1 - mW3 , 

m=0 n=0 

oo m 

^ = EE^A 4(1 - ro) " /3 , 

m=l n=0 

oo m 

m=l n=0 



(3.12) 



This particular ansatz will not be sufficient to include all cases; for example we have seen in 
section 3.1 that the GTV solutions with h\ = 2N C have a completely different form (3.5) in 



the UV. We will find other limits in which this is the case, and which will have to be treated 
separately. 

As in the IR, we substitute the ansatz (3.12) into the equations of motion, and in this 
case we find nine independent parameters, which we take to be 

K 00 , K 30 , H 10 , H n , $io, $30, W 20 , W m , U40, (3.13) 
and which we again relabel to make contact with the SUSY case: 



00 



2c+ _Qo _ 



64e 4 Poc 3 



30 



48c 2 + 

The nine relabeled independent parameters are then 

c+, c_, $oo, Q , p , H u , W 2 o, $30 
and the expansions are 

2g _ „ Ap tAU „ , ^ , o„ w2 \ , /n/'„-4p\ „2/i _ C + j 



W 4 o = 2e po . (3.14) 



U 40 , (3.15) 



„2ft 



(4H 11 p + Qo + 2c+Wi ) + O(e—3f), e 



3 3 

_ 2 

a = W2oe sP + 
, 9TU 20 



2c\ r ~ I 
e -2p + ( e -fp) 5 



Qo 
4 



4<£>oo $30 



-e 3^ + 



+ f 4e^ - _ 2W%\ + 1 g _ 2p + 

3 V c + 6 / J 



(3.16) 



e 3 



The most significant difference here when compared to the SUSY expansions (2.15) is the 
presence of the new terms at leading order in the UV in a and b. This corresponds to the 
presence of the additional terms proportional to p -1 / 2 in a which we saw in the GTV solutions 
(3.5), and we will see in section [5] that the interpretation as a mass term still applies. The fact 
that the extra terms we obtain here are exponential rather than polynomial in p is related to 
the qualitatively different UV asymptotics in the baryonic branch (2.15) as opposed to the 
CVMN solution (El9b. 
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We can recover the SUSY case from (3.16) by setting 

N c 

Hu = ^, W 20 = 0, $ 



3N ( 



30 



4rf 



c e 4* c 



(3iV c + 4Q ), V, 



40 



e 2p °(N c + Q ). (3.17) 



For the regular SUSY solution (2.14, 2.15) we also need p Q = and Q Q = —N c . 

In summary, our solutions are described by the fourteen parameters: the five from the IR 



(3.10) and nine from the UV (3.15). However, if we consider only solutions which match both 
the IR and UV expansions (3.8, 3.12) these are clearly not all independent. There can be at 
most five independent parameters, as the required solutions can be parametrised by the IR 
boundary conditions alone. However, we generically expect even fewer. 

Our goal is to find a solution which smoothly interpolates between the IR and UV expan- 
sions. This will require that these two parametrisations lead to identical functions. We can 
express this as a system of twelve equations Q 

g{hi ...w 2 ;p) = g(c + . . .V i0 ;p), ^d(h ■ ■■w 2 ;p) 
h(hi ...w 2 ;p) = h(c + . . . V A0 ;p), ^K h l ■ ■■ w 2',p) 



d_ 

dp 



g(c+ . . .V A0 ;p), 



d h(c + ...V 40 ;p), 



dp 



(3.18) 



b(hi ...w 2 ;p) = b(c + . . . V40; p), 



w 2 ;p) = f- p b(c + ...V m ;p). 



This system can be further reduced using the constraint (appendix[A]). This means we can for 
instance express the derivative of one of the functions in terms of the other functions and their 
derivatives. This leaves us with a system of eleven independent equations which we would 
expect to allow us to solve for eleven of our fourteen parameters. Although in principle further 
redundancy in the system of equations (3.18) would allow for more independent parameters 



up to a maximum of five, the numerical analysis discussed in |18j and below appears to 
support this conclusion. Of the three remaining parameters, one corresponds to our ability 
to shift the dilaton, which has no other effect on the solution. The final two parameters we 
then associate with movement along the baryonic branch and finally the breaking of SUSY. 

In much of the following it will be convenient to describe the solution space in terms of 
the parameters that appear in the IR expansions. Firstly, as in [18] the smaller number of 
parameters makes finding suitable numerical solutions much simpler starting from the IR. 
Secondly, our IR ansatz (3.8) imposes a comparatively natural restriction on the solutions, 
while the UV ansatz (3.12) is more arbitrary, merely being a plausible candidate for a gen- 
eralisation of the most usual SUSY solution. Indeed, as discussed we know that it does not 
apply in several interesting special cases. 

To allow contact with the SUSY case, we choose h\ to parametrise the position along the 
baryonic branch. We could then in principle choose any combination of the remaining IR 
parameters v 2 , w 2 and k 2 to describe the remaining degree of freedom (figure [6]) . It turns out 
that a description in terms of v 2 is usually simplest; we see from (3.9) that its SUSY value, 
-2/3, is independent of h%. 



v: 



SUSY 



7 We write the functions resulting from a given choice of the IR parameters {hi, v%, W2} in the form 
g(hi,k2,V2, W2',p). Similarly the expressions of the form g(c + , c_ , Q , po, Hu, W20, $30, V40; p) refer to the 
functions resulting from a given choice of the UV parameters. 
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Figure 6: The space of solutions, seen in terms of (v2,h±), (w2,h±), and (&2, /ii), where 
^2 = kijhx. The blue curves denote the SUSY baryonic branch (section , and the green 



lines correspond to the GTV solutions (section 3.1). Note that all the GTV solutions have 
&2 = 0, as can be seen from (3.3). Any solution to the equations of motion of the form (3.8) is 
represented by a point on each of these diagrams. If we require a well-behaved UV, specifying 
the position on one diagram is sufficient to determine the positions on the other two. For 
example the marked point represents schematically a generic solution of the sort presented in 
|18j . The values marked at the top show the SUSY values in the limit h\ — >■ 00, corresponding 
to the KS solution [2]. 



3.3 Finding globally regular solutions 



In order for us to be able to conclude that the IR expansions of the form (3.8) and the UV 
expansions of the form (3.12) describe the same system of solutions, it is necessary to find 



numerical solutions interpolating between them. This was achieved in |18j for isolated exam- 
ples, simply by manually searching the IR parameter space for solutions with the expected UV 
behaviour. However, without having a good understanding of the structure of the parameter 
space it was difficult to make progress. In particular, the approach was in practice limited to 
solutions very close to the SUSY case (i.e. V2 ~ —2/3). 

Fortunately, we can make use of the simpler system of GTV solutions (section 3.1). Just 



as the CVMN solution (2.17) can be obtained from the SUSY baryonic branch solution (2.15 
2.14) in the limit hi — > 2N C , we would expect to obtain the GTV solutions from our non-SUSY 



generalisation of the baryonic branch in the same limit. In the IR, this is indeed the case; by 
setting W2 = V2 and &2 = in our solution we recover (3.4). Of course, there is no way to 
obtain the GTV UV (3.5) from UV expansions of the form (3.12), but this is to be expected 
given that the equivalent statement is also true in the SUSY case — the CVMN UV (2.19) 
cannot be obtained as a simple limit of the generic UV (2.15). 

As the GTV system has no redundant parameters in the IR, it is simple to generate 
numerical solutions. It is then possible to deform this well-understood case by increasing h\ 
and adjusting W2 and V2 slightly to correct the UV behaviour. More precisely, for a given 
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value of V2 = Av2 — 2/3 it is trivial to obtain a numerical solution with hi = 2N C , for which 
W2 = V2 and = 0. We then deform this by keeping Av2 fixed and setting hi = 2N C + Ahi. 
If we use a small perturbation Ahi, we will require corrections of the form 



W2 = w| UbY (A/ii) + Av 2 + 6w 2 (Ahi,Av 2 ), k 2 = ^ {Ahi) + 5k 2 (Ah 1 , Av 2 ), (3.19) 
where 5w2 and <5fc2 are extremely small. 



The far UV of the solutions obtained in this way match our general ansatz (3.12 ), justifying 
our assumption that the GTV solutions can be viewed as a limit of our deformations of the 
general case. 

In itself, this yields a considerable advance over using only the approach described in 
[T%] - - it gives us access to solutions with hi fa 2N C and general v 2 , in addition to those 
with V2 ~ —2/3 and general hi. More significantly, however, it allows us to understand the 
behaviour of solutions with generic values of both hi and V2 in terms of the corresponding 
solutions in the two limits. 



4 The two-dimensional solution space 
4.1 Combining the effects of hi and v 2 



As we have seen in section 3.2 , the system is described by a two-dimensional parameter space, 



corresponding to the position along the baryonic branch and the size of the SUSY-breaking 
deformation. We generate numerical solutions starting from the IR, so we are led to the 
choice of hi and one of the three SUSY-breaking parameters {u>2, V2, A^}. Of these V2 turns 
out to be most convenient because u| USY is independent of hi . 



In section 2.3 we described the effect of varying hi in terms of the scale correspond- 



ing to the transition between the CVMN behaviour (2.17-2.19) and the generic (KS-like) 



behaviour (2.15). Similarly, in section 3.1 we introduced the scale psvsY: associated with the 



transition between the qualitatively SUSY CVMN behaviour and the (non-SUSY) GTV UV 



(3.5) 



In the case of a generic solution, with hi > 2N C and i>2 ^ —2/3, we find that these features 
survive and both scales are present. The sequence then depends on the ordering of the two 
scales. If ph x < psiJSYj the sequence is (figure [7] (i)): 

-2p 
-2p 
-2p/3 



P < Phi 

Phi < P< PSUSY 
P > PSUSY 



~ 9 ~ 
~ 9 ~ 



constant, 
2p/3, 
2p/3, 



a fa b 
a ~ b 
a ~ b ~ e 



~ e 
~ e~ 



(SUSY, CVMN-like) 
(SUSY, KS-like) 
(non-SUSY, KS-like) 



On the other hand, if psusy < Phi we have (figure [7] (ii)) 



P < PSUSY 
PSUSY < P < Phi 

p > phi 



a 

a m 
a ~ 



-2p 



k ~ g ~ constant, 
k ~ g ~ constant, 
k ~ g ~ 2/9/3, 

It appears that ph x is almost independent of V2, and that psusy is almost independent of 
hi, although this may break down for sufficiently large hi and t>2, depending on the precise 



b ~ e 
b ~ p- 1 / 2 
b ~ e- 2 ^ 3 



(SUSY, CVMN-like) 
(non-SUSY, GTV-like) 
(non-SUSY, KS-like) 
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definition used for the scales. In fact the presence of the two scales becomes less clear as they 
move into the IR for large hi and v 2 . This reflects the reduced gradient far from the CVMN 
solution {h\ = 2N C ,V2 = —2/3) in figures [2] and [4j We show the behaviour of the functions 
for some generic solutions in figure [8} In this case h± is large enough that ph x is not visible. 




Figure 7: Plots of log a against p comparing the solutions obtained for each combination of 
hi = 2N C , hi = 2N C + Ahi, v 2 = -2/3 and v 2 = -2/3 + Av 2 . 

(i) Ahi = 1CT 5 , Av 2 = 1CT 9 , with N c = 1. Here p hl < Psusy- In the IR we see the CVMN- 
like behaviour, with a ~ 2p/ sinh2/3. At ph x the solutions with hi > 2 deviate from this, but 
after the transition the gradient is unchanged as we still have a ~ e~ 2p . Then at psuSY the 
non-SUSY solutions switch to the slower decaying behaviour. 

(ii) Ahi = 10~ , Av 2 = 10~ 4 , again with iV c = 1. Here psusy < Phi- The IR still shows the 
CVMN-like behaviour. At psusy the non-SUSY solutions switch to the GTV-like behaviour, 
with a ~ p~ l l 2 . Then at ph ± the solutions with hi ^ 2 show a transition. In the SUSY 
case the gradient is the same after the transition, but in the non-SUSY solution the gradient 
increases. This corresponds to the transition between a ~ p -1 / 2 and a ~ e~ 2p / 3 . 



4.2 The boundaries of the parameter space 

A notable feature of the GTV solutions is the restriction to —2 < v 2 < for solutions with a 
regular UV. There is no obvious way to determine whether an equivalent condition holds for 
hi > 2N C , or to find the correct generalisation. 

However, numerical observations suggest that w 2 (hi, v 2 ) becomes independent of hi for 
for v 2 — > 0, and that in particular there is a family of solutions with a = b = 1 and g = k 



even for hi > 2N C . This corresponds in our IR expansions (3.11) to setting 



hi AN 2 5,c TI c Y fl . 

w 2 = v 2 = 0, k 2 = ^-—± = -A| USY , 4.1 

3 Shi 6 

which agrees with the values obtained numerically. Setting v 2 > (so that a > 1 for small p) 
appears numerically to result in a divergent UV, and it seems likely that this is indeed the 



correct generalisation of the boundary. This corresponds to the solid red curves in figure 10 
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(i) (ii) 

k k' 




P 




Figure 8: Plots of some of the metric functions for different V2, having set h\ = 2.3 with 
N c = 1 and 0o = 0. In (i) we plot k, showing that the SUSY-breaking parameter V2 has little 
effect on the qualitative behaviour except in the case vi = 0. This transition between the 
generic k ~ 2/9/3 and k ~ p/V% is shown clearly in (ii), in which we plot the derivative. The 
UV behaviour, and effect of V2, is very similar in g and h. In (iii) we plot the dilaton, showing 
that is a function of V2 for constant hi, and (iv) shows the effect on a. The values of V2 
used, and the corresponding values of W2 and k<z, are shown as coloured points on the solution 
space diagrams in (v). The colours correspond to those of the curves in (i)-(iv). As in figure 
6j in (v) the blue curves are the SUSY solutions and the green lines are the GTV solutions. 
The shaded areas are the regions which are excluded according to the discussion of sections 
4.2 4.3 As we will explain in section 4.3, we can restrict our attention to W2 > —2 without 



loss of generality. 



18 



Setting a = b = 1 and g = k in the equations of motion, we find that our UV ansatz 
(3.12), (expansions in powers of e 4/9 / 3 ) is not suitable. However, using equivalent expansions 



in powers of e^ p does lead to a solution: 

'^00-^20 + ^, 



,:2ii 



2k 



2 + 



K 00 e^ zp +{K 20 + 



2K 00 

iV 2 



+ 



oo 



2V2K o P 

e -V2 P + ( e 



-V2p 



+ 0(e- 3 ^) 



1 



1 



K 2 



4K 00 ^20 + iV c 2 + 2V2iV c 2 p 



2 ^ + Ofe 



(4.2) 



It is important to emphasise that although the form of (4.2 ) is simply the original ansatz (3.12 ) 
with the replacement 4p/3 — > \/2p, we cannot obtain these expansions from the generic UV 
(3.16) simply by a change of coordinates. For example, here we have e 2k = e 2g , whereas in 
(3.16) we have e 2k ~ 2e 2fl /3 for large p. This is why we have not attempted to match the 



parameters in (4.2) to the usual set {c+, c_, . . . }. Instead, we denote the two free parameters 
by Kqq and K20, the leading parameter (roughly corresponding to c+) being Kqq. Note that, 
as we have set v 2 = 0, the two parameters Kqq and K 2 q cannot be independent once we match 
to the IR. This is analogous to the SUSY solutions, in which there are two UV parameters c+ 
and c_, which are related by the requirement to match to the (one-parameter) IR solutions. 



In section |3.1 1 we noted that the 'twist' which mixes the S 2 and the S 3 could be removed 
by a change of coordinates when a = b = 1. As we still have g = k here, the same coordinate 
transformation still works, leading to a simplified system. With C\ and F5 unchanged from 



(2.9), we now find 



ds% = < 



0/2 



■2k 



h^dx 2 ^ + h 1 ' 2 ( e 2k dp 2 + e 2h dn 2 + ^dfi 3 



N, 



c Cj\ A Cj 2 A Q3, 



H 3 = 2Ar ce 2?i - 2fc + 2 *-*- s[q9 dp A dQ A d(p _ 



(4.3) 



Unfortunately the boundary for v 2 < —2/3, corresponding to v 2 = — 2 in the GTV solu- 
tions, seems to be much less accessible numerically, in part due to the presence of changes of 
the sign of a and b. However, in the next section we will shed some light on this matter. 



4.3 A Z 2 symmetry 



The system we describe in section 2.1 which applies to all the solutions we consider, exhibits 
a 7L 2 symmetry X which exchanges the two 2-spheres of the conifold and changes the sign of 



the 3- forms F3 and H$ in (2.9). To see this, we make use of the fact that all the systems we 



consider can be described by the Papadopoulos-Tseytlin ansatz [36]. This can be written in 



19 



the form 



ds 2 E = e^ 2 (h~ 1 / 2 dx 2 3 + h 1 / 2 ds 2 ), 



ds P 



-Sp+3q 



(4dp 2 + g 2 ) + e 2p+3q \ cosh 



where the angular forms g§ and Ui are given by 
U\ = d9, UJ2 = — sinQdip, 



e z {<jj\ + uj 2 ) + e z (Cj\+CjI 
— 2sinhy (oj\Cj\ + 002^2 

55 = W3 + cos 9 dip. 



(4.4) 



(4.5) 



We use here the notation of |17j . in anticipation of making contact with their results in 
section 4.4| 8 By comparing the metrics in the two cases we can write an explicit relation 



between our original functions and those used in (4.4): 
4 



,10p 



a g+h-2k 



Ag+Ah+2k 



-h 



V / 4e 2/l 



e 2 9a 2 



ae- 



e -9y/4e 2h + e 2 9a 2 . 



(4.6) 



It is then possible to show that the metric and fields are unchanged (up to a change of 
sign) if we exchange (8, <p) o (6, <p) and relabel z o —z. In the KS solution [2] which we 
obtain by taking the limit h\,c + — > 00 in the SUSY solutions (section [2]) , z = and the 
transformation X reduces to a simple change of coordinates. This is the iVf = version of the 
Seiberg duality discussed in [281 EH ED] • 

We now consider the effect of X on a generic globally regular solution of the sort we have 
discussed, for which z 7^ 0. Inverting (4.6) we find 

e 2g = 2 6 / 5 e 2p+3q - z coshy, 
2 ll/5 



e 2fc = 2 -4/5 e 2 P +3 g+ * sech 



V, 



/2k 



-8p+3q 



It is then clear that the effect of taking z 



a = e tanhy. 
-z can be written as 



,2 S 



^g+2z 



J2h-2z 



a 



e- 2z a. 



Referring to our expansions (3.11 3.16), we find 

1 + (2 + w 2 )p 2 + 0{p A ) 
1 



for p 



1 + — ^4ffnp + Qo + ^c + W 2 z ) e~ 4 "/ 3 + 0(e- 4 "/ 3 ) for /» 





00. 



(4.7) 
(4.8) 

(4.9) 



This means that the transformation (4.8) has only subleading effects on g, h and a, and in 
particular the transformed functions are still compatible with the form of our expansions 
(3.8, 3.12). More specifically, we can see from (4.9) that for z —> —z we need to take 



Wo 



w 2 , 



(4.10) 



We adapt the notation slightly to avoid confusion with the vielbeins (|2.8||. The relationship with [T7] is 



.here 1171 „ ] ~here 

Ji — e ■ and w,- 



We also have phere = TTf7]/2. 
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2. As implies that w| USY 



corresponding to a reflection in the line W2 - 
h± — > oo, this is compatible with the fact that the KS solution has z 



-2 for 



0. Note that because 



k and b do not appear in (4.8) we can conclude that /c2 and V2 are unchanged under X. 



This gives us a simple procedure whereby for each solution discussed in 4.1 specified by 
values of (hi, V2), we can obtain different solution, with a different value of u>2- Because of the 
way that I acts only on the subleading terms in the UV expansion we can be sure that the 



'reflected' solution will also be globally regular and compatible with our ansatz (3.12). This 



can be seen numerically. For a given pair of values of (hi,V2), the requirement of UV regularity 
gives us values of W2 and &2 as described in section |3.3| If we then take W2 — > —4 — W2 we 
immediately find another solution with the correct UV behaviour, without having to adjust 
t>2 or &2- The two types of solutions are compared in figure [9) 



(1) 




di) 




(IV) 




Figure 9: Comparison of the solutions before and after the transformation (4.8). The blue 
solid curves correspond to the original description with W2 > 
to the 'reflected' solutions with W2 < 



-2, and the red dashed curves 
-2. In (i)-(iii) we plot the three functions {g,h,a} 
which are affected by the transformation, and in (iv) we show z, as defined in ( |4.6| ), for which 
the transformation is simply a change of sign. These solutions have hi = 2.3 (with N c = 1) 
and V2 = —1/3, resulting in &2 ~ 0.195 and W2 ~ —2 ± 1.58 (corresponding to the yellow 
plots in figure [8]) . 



We should emphasise, however, that although this results in a distinct solution to the 



equations of motion (A.2-A.7), it does not actually correspond to a different background 



X is simply a relabeling, which is obscured by our choice of basis for the functions. 
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By demanding that z 



-z while the other functions are unchanged, we can also write 



down the effect on the UV parameters equivalent to (4.10): 



Qo 
c+ 



:»2U 20 ( 2— +3W 2 2 



Qo 
c+ 



3Wi , 



2e 



2po 



Qc 



w 20 



-w 



20 



Qo 



W 2Q 



H 



11 



c+ ' 



(4.11) 



while keeping the remaining parameters {c+, $oo, W20, V40} fixed. We retain the factors 
of l/c+ in anticipation of taking the limit c+ — > 00. 

For solutions with W2 = —2, the IR expansion (4.9) appears to vanish at all orders. We 
would therefore expect that these solutions have z = for all p, meaning that as in the KS 
solution X is a symmetry of the geometry. Our numerical calculations support this assumption 
— for these solutions we find that z is indeed essentially zero everywhere (we find z < 10~ 
for all p < 30). 



This family of solutions consists of a line in the (hi, W2) plane (see figure 10), and it would 
be interesting to determine the corresponding curves in the (/ii,t>2) and (^1,^2) planes. We 
have not been able to determine exact expressions for these functions, but for large hi we 
find numerically that A-t^/ii) = ^2(^1) + 2/3 ~ — 1/hf, and 



Ak 2 (hi) = k 2 (hi 



fc 



SUSY 



(h 



16 
45/ii 



e(hi), 



(4.12) 



where the higher-order corrections e > to this last expression are extremely suppressed. 
For example, with hi ~ 10 3 we find that using Ak2 = 16/45/ti gives the correct value up to 
around eleven significant digits. In figure 10 the curves ^2(102 = —2) and ^2(^2 = —2) were 
obtained from expansions in powers of I//11 fitted to eight solutions determined numerically. 

If these solutions are indeed symmetric under I for all p then we can write down a 
relationship between some of the UV parameters, analogous to the requirement that W2 = —2. 



Specifically, referring to (4.11), we fin 



W. 



20 



2Qo 

3c+ ' 



H 



11 



0. 



As expected, this is satisfied by the SUSY values (3.17) in the limit c + 
to the KS solution. 



(4.13) 



00, corresponding 



In section 4.2 we considered the generalisation to hi > 2N C of the upper bound v% = 



in the GTV solutions (section 3.1). It is suggestive that the line of solutions in which the 
geometry possesses a Z2 symmetry passes through the lower bound, V2 = —2. In the light of 
the discussion in this section, we should reinterpret this boundary in the GTV solutions. If 
we parametrise the solutions by w%, we see that there is no lower bound on W2, but ^2(^2) 
has a minimum at W2 = —2. This description was not possible in the context of [9], in which 
all solutions had a = b (so that V2 = 102)- 



9 Of course, we still have the usual undetermined relationships between the UV parameters, so that we are 
left with only one degree of freedom corresponding to the position on the line W2 = —2. 
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Interpreting the boundary as a minimum of ^2(^2) would imply that the line W2 = —2 is 
the right generalisation to h\ > 2N C . This is supported by our numerical analysis. It appears 
not to be possible to tune to a regular UV for values of V2 smaller than that which gives 
w 2 = -2. 

Of course, we must be cautious here — our inability to find a solution with u>2 < — 2 could 
simply be the result of a significant discontinuity in the values of the other parameters across 
the line W2 = —2. 





1 10 F 




q 10 F 



-4 -10/3 -8/3 -2 -4/3 -2/3 

w 2 




1/10 1/5 3/10 
k 1 =k 2 /h l 



Figure 10: The space of solutions, as in figure [6j Again, the blue and green curves are 
the SUSY and GTV solutions respectively. The red curves correspond to the case a = b = 



1 discussed in section 4.2, while the orange curves correspond to the solutions which are 



invariant under I, and so have a Z2 symmetry of the geometry (section 4.3). The dotted 
curves are the equivalents with W2 — > —4 — W2- Under the assumption that these two cases 
constitute the correct generalisation of the requirement —2 < V2 < in the GTV solutions, 
the gray shaded areas show the regions where no regular solutions exist. 



4.4 The limit h\,c + — > 00 

Having discussed a non-SUSY generalisation of the baryonic branch, it is natural to consider 
the generalisation of the Klebanov-Strassler solution [2j itself, which in the SUSY case occurs 
in the limit h\ ~ c+ — > 00. In terms of the functions {p, q, y, z} which we introduced in section 



4.3, the SUSY KS solution has a simple exact description: with <I> = constant and z = 0, we 



have 



o5/4 

e i0p = R 3 sinh2 p, e 15q = -^K 2 smh i 2p, e y = tanhp, (4.14) 



where we have defined 



(smh4p-4p)V3 

2 1 /3 s i n h2 / o V ; 
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The remaining function b = 2p/ sinh2 / o is the same as in the whole SUSY baryonic branch. 
As we have seen in section |4~3| the fact that z = implies that the geometry itself possesses 
a Z2 symmetry. 

Of course, in order for the concept of a non-SUSY generalisation to be meaningful, we 
have to choose which characteristics of the SUSY KS solution we want to keep in the non- 
SUSY solution. One natural possibility would be to require that the geometry retains the 
Z2 symmetry, in which case the we obtain the family of solutions with w 2 = —2 which we 
discussed in section 14.31 

However in [T7], Dymarsky and Kuperstein (DK) followed a different approach. They 
noted that the KS background has several simplifying features which are retained in the 
linear deformations studied in |3T|, [38], but not in the generic baryonic branch: 

(i) A constant dilaton, e* = g s 

(ii) An imaginary self-dual 3-form fluxp'j 1G3 = *6G3, where G3 = F3 + jj-Hz 
(hi) An RR 4-form satisfying C 4 = if^Voli^, where ds 2 = H~ l / 2 dx 2 ^ + H^ds^ 
(iv) A Ricci-flat 6d unwarped metric 

As noted in [IT], these are particularly convenient because they mean that the fluxes com- 
pletely decouple from the equations which determine the metric. It should be noted that in 
our solutions (ii) and (hi) are satisfied automatically once (i) is imposed. 

By imposing that these properties are retained, DK found a one-dimensional family of 
solutions which break both SUSY and the Z2 symmetry of the geometry (although the full 
symmetry including the exchange z -H- —z is of course retained). It seems natural to assume 
that this corresponds to a line of solutions in the two-dimensional solution space described 
above. 

To see that this is indeed the case, we first need to identify the appropriate limit. Referring 



to our generic IR expansions (3.11), we see that we obtain a constant dilaton in the limit 
h\ —7- 00, as in the SUSY case. This means that conditions (i)-(iii) are satisfied. It is also 
possible to check that this results in the IR expansion for the 6d Ricci scalar vanishing, as 
required by condition (iv). 

We now look to relate our three SUSY-breaking parameters {w 2 , k 2 ,v 2 } to the parameters 
{Ci> C25 C3} used in [T7]. Looking then at the IR expansion for z, we find by substituting our 



IR expansions (3.11) into (4.6) 



z = {2 + w 2 )p 2 + 0(A (4.16) 
meaning we can compare with the expression given in [T7] and conclude that 

w 2 = 4Ci - 2. (4.17) 



3 Here *6 is the six-dimensional Hodge dual 
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To gain the relation for fc 2 we look at the expansion for e y and upon taking the limit h\ — v oo 
we find 

e» = p-Q+4C?)p 3 + 0(p 5 ). (4.18) 

This does not have enough freedom in the p 3 term when compared to |17j . To fix this, it is 
possible to take fc 2 — > oo while keeping fixed &; 2 = h%lh\. This then gives 

^+4C 1 2 -^ 2 )p 3 + 0(p 5 ), (4.19) 

which we can match to the result of [T7] by setting 

^2 = ^ = ^(13-90C 2 ). (4.20) 
h\ bo 

We finally need to determine the relationship between V2 and (3. This can be achieved by 
comparing our expansion for b with that for F = (1 — b)/2 in [T7], from which we obtain 

«2 = -|(C3 + l). (4-21) 

In summary, in the limit h\ — > oo we find the following relationships between our three 
SUSY-breaking IR parameters and those used in [T7] : 

w 2 = Ki- 2, k 2 = ~ = ~(13-90C 2 ), «2 = -?(C 3 + l). (4.22) 

All DO 3 



Of course, setting the £i to zero we recover (the large-/ii limit of) the SUSY values ( |3.9[ ). 
In fact, defining for example Aw 2 (hi) = W2 — w 2 (a*i), we obtain 

1 13- 3 

Ci = -A™ 2 , C2 = --AA: 2 , C3 = -7 > A^ 2 . (4.23) 
4 3b 2 

In the UV we are less sure how to find similar relationships between parameters. It is 
clear from the numerical analysis that the relevant limit is still c+ — > oo (even if the precise 
relation (2.16) may no longer hold in the non-SUSY case), and we know we will need $oo —> (f>o 
in order to get a constant dilaton. However, it is not obvious how the other parameters in 
(3.15) behave in this limit. One possibility is suggested by the fact that in the case of the 
IR parameters we could have guessed the correct behaviour (u 2 ~ io 2 ~ constant, k 2 ~ hi) 
from the h\ -dependence of the SUSY values (3.9) in the limit. Using the same approach in 
the UV would imply that we should consider all the remaining parameters fixed except for 
$30 ~ l/c+. 

Looking at the UV expansions for the 6d Ricci scalar and the dilaton we find that in fact 
the limit c+ — > oo is itself sufficient for Ricci-flatness, and taking both c+ — > oo and —> 
gives a constant dilaton. This can be seen for the SUSY baryonic branch in figure [l] (iv); the 
non-SUSY solutions show qualitatively the same behaviour. 

Unlike in the case hi —> 2N C , our numerical approach does not allow us to take the limit 
hi — > oo explicitly. However, we can probe sufficiently large values of hi to yield solutions 
which appear to have many of the characteristics we expect from the true limit. For example, 
we do not have to take hi very large before the dilaton is very close to constant. Notice in 
figure [T] (iv) the curve for hi = 12 appears to lie on the axis. 
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5 Remarks on the dual field theory 



Here we shall discuss a little about the dual field theories to the gravity backgrounds we have 
presented. Much of the following is similar to that of |18j . although we have not restricted 
ourselves to small deformations, as was discussed in that paper. 

We will only consider the solutions with hi > 2N C . In this case the geometry is 'almost' 
asymptotically AdS§. More precisely, for large p we can write the metric in the form 



ds 2 ~ 



j2 H{u) 1 ' 2 j 2 j2 

"> x l 3 H o +« s 5i 



H(u) 



H{ u yi 2 -~ it- 
log u + constant + 0(u~ 2 



(5.1) 



where we have defined a suitable radial coordinate u (increasing with p). For the generic 



solutions satisfying the ansatz (3.12) (including the SUSY solutions), the definition which 
results in (5.1) is u = e 2p / 3 . For the solutions with v 2 = discussed in section 4.2 we instead 
need u = e p l^ 2 . The term of order log-u in the correction H(u) results from the sub-leading 



behaviour of the dilaton (3.16, 4.2). 



There are three different field combinations which are invariant under the rotation which 
are of interest 1211 . The first is the dilaton and the others are defined as 



M 



„2z 



1 = a 2 + 4e 2h - 2 ° - 1, 



Mo 



2h+2g-4k 



(5.2) 



In the case of the generic solutions described by (3.16), these functions have UV expansions 

'3N 2 



1 



2ct 



-P 



-4$ c 



$30 

4 



Mi = (SH llP + 3c + Wi + 2Q ) 
9 27 



e -fp + ( e -4p 

,-4p/3 



Mo 



:W 2 e 



-4p/3 + 0(e -8p/ 3)> 



+ 0(e 



8 P /3n 



16 16 1 ' (5 ' 3) 

By looking at the asymptotic behaviour of fields (and combinations of them) it is possible 
to think of our constants in terms of the operators which are deforming a fixed point. We may 
do this as it is understood that a generic field A4 ~ as u — > oo behaves in the following 
manner. If A > (or A = 0) it is either an indication of a relevant (or marginal) operator 
in the Lagrangian or the VEV for an operator of dimension A. If instead, A < 0, then it 
indicates the insertion of an irrelevant operator of dimension (4 — A) in the Lagrangian. 

Using this analysis it can be seen, from the UV expansion above, that the dilaton falls 
into the marginal operator category as it has scaling dimension A = 4 (this can be associated 
with a certain combination of gauge couplings discussed in |18j). 

We can further use this analysis on the expansion of the function b(p) presented here for 
convenience 



91U20 



4p 



+ 



lOWfn o 



4e 2p ° 



Q0W20 



23Wf 



p + V40 



e~ 2p + 0(e~ir p ). (5.4) 
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Here we can see that W20, which we could consider to be our 'SUSY- breaking constant', 
corresponds to an operator of dimension three being inserted in the Lagrangian. We antici- 



pated in section 3.2 that we can associate this operator with the mass of the gaugino, as in 



]. Following the SUSY case we also associate e 2po , which appears at next-to-leading order 
in M\, with the VEV of the gaugino. From this we can write schematically 



W20 — > m\\, 



= 2p 



(AA) 



A 3 

A YM- 



(5.5) 



It should be noted that this association is not exact once we have broken SUSY — the 
SUSY-breaking parameter can generically also deform the gaugino VEV, as indicated by the 



contributions from W20 and V40 to Mi in (5.3). 



As discussed in appendix [Bj it appears that W20 — > 00 as we approach the boundary at 
V2 = W2 = (figure 11). This suggests that we can interpret the solution on the boundary, 
with a = b = 1 for all p (section 4.2) as corresponding to a field theory in which the gaugino 
has been given infinite mass. We therefore no longer have soft SUSY breaking — the theory is 
non-SUSY all the way into the UV. Presumably, by sending the mass to infinity we effectively 
remove the gaugino entirely, obtaining a completely non-SUSY theory. 

We can now look to the field combination Mi and see that it can be thought of as corre- 
sponding to the VEV of a dimension two operator U. In the SUSY case we can identify [39 j 



Z^~tr [AA*-B*B], 



(5.6) 



and this operator getting a VEV is the exact thing which allows us to explore the baryonic 
branch. Notice that in the SUSY case W20 = and the leading term of Mi vanishes for 
c+ — > 00, when we recover KS. This is also the limit in which the geometry is invariant under 



the Z2 symmetry X which we discussed in section 4.3 In fact, from the point of view of the 
field theory, the transformation X can be identified with swapping A B |37j . 

As soon as we move away from the SUSY solutions we can no longer make the identification 



(5.6). However, it is still instructive to consider the behaviour of the operator IA associated 



with Mi. From (5.3) it is clear that we can expect IA to be changed when we break SUSY 
while keeping c_|_ fixed. Indeed, referring to the definition (5.2), we see that Mi = when 
z = 0. This applies at all p in all the solutions on the line W2 = — 2. (As required, we 



see that the combination of parameters appearing in the UV expansion (5.3) vanishes when 



(4.13) is satisfied.) It is interesting that the presence of the Z2 symmetry still corresponds to 



the vanishing of this operator, even in the non-SUSY case. This is perhaps indicative of the 
extent to which the structure of the SUSY system survives in the generic case. 

As we move in the opposite direction from the SUSY solutions, increasing V2 (and W20) 
we find numerically (appendix [B]) that both terms at leading order in Mi diverge. However, 



in the limit we obtain the solutions described in section 4.2 and the expansions (5.3) are no 
longer valid. Instead, for large p 



Mi 



2iV c 2 
^00 



+ 0(e 



(5.7) 
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This is qualitatively different to the generic case. Firstly, we now have M\ — > 2 in the UV, 
as opposed to Mi — > 0. This indicates that these solutions do not recover the Z2 symmetry 
in the UV. Secondly, the next-to-leading term is now of order u~ 4 , meaning that we can no 
longer associate this field with a dimension two operator. 

There is some subtlety here in the fact that unlike in [18] we have allowed our deforma- 
tions of the SUSY solutions to become large. It is then not clear that any deductions based 
on analogy with the SUSY solutions remains valid. In particular, we cannot not necessarily 
expect to find stable solutions for all values of W2o- However, the similarities between the 
SUSY and non-SUSY solutions are interesting. It should be noted that we still find a con- 
tinuous and smooth deformation of the SUSY solutions between smaller and larger values of 
the non-SUSY deformations in the IR. We only find a different UV expansion in the limiting 
cases (or boundaries of our solution space). 



6 Summary and conclusions 

In this paper we study the full two-dimensional space of solutions which can be considered 
to be the non-SUSY generalisation of the baryonic branch (extending the work of [H]). We 
include the solutions compatible with the PT ansatz which have both a regular IR, of the same 
form as that of the baryonic branch, and are related to the baryonic branch by a continuous 
change of parameters. 

In addition to the SUSY baryonic branch and its limiting cases (Klebanov-Strassler and 
Chamseddine-Volkov/Maldacena-Nuhez), this solution space also includes two previously 
studied one-dimensional families of non-SUSY solutions as limits. In the limit which yields in 
the SUSY case the CVMN solution we obtain the solutions of Gubser, Tsyetlin and Volkov 



[9] (presented here in section 3.1), while in the limit corresponding to the KS solution itself 



we obtain those of Dymarsky and Kuperstein [T7j (presented in section 4.4). The behaviour 
of generic non-SUSY solutions lying away from these boundaries can be understood as a 
combination of the effects which are present in the SUSY baryonic branch and the GTV 
solutions. 

Alongside these cases we identify two additional one-dimensional families which are of 
interest. The first is the boundary of the solution space with v 2 = W2 = corresponding 
to the positive boundary of the GTV solutions. Here we can no longer argue that SUSY is 
softly broken (the gaugino mass appears to be infinite), and we find that a = b = 1 for all 
p. Notably, this changes the geometry to an explicity non-SUSY case (a cone over S 2 x S 3 ). 
We also find an explicit UV expansion for the solutions on this boundary which is different 
from the generic UV. The second family lies on the line W2 = —2, upon which the geometry 
possesses a Z2 symmetry just as in the Klebanov-Strassler solution. This family of solutions 
corresponds to the other boundary of the GTV solutions. 

Moving away from the boundaries, we have also shown that solutions with W2 < —2 are 
related to those with W2 > — 2 by a Z2 symmetry and describe the same physical system, 
although the solutions themselves appear different. In the two-dimensional solution space 
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much of the SUSY structure survives. In addition to the various quantities calculated in [18J, 
which are mostly unaffected by SUSY-breaking at leading order, we find that the presence of 
a Z2 symmetry of the geometry is still linked to the vanishing of a dimension-two operator. 
In the SUSY case this reflects the fact that in the dual field theory the Z2 transformation 
corresponds to the ability to interchange the baryons. 

It would be interesting to know to what extent this description applies to the non-SUSY 
case. To address this, it would be necessary to gain a more detailed understanding of the 
field theory, including calculation of the mass spectrum. Another question which we did 
not address is the issue of stability. It would be useful to determine if, and how much, 
the parameter space is restricted by this requirement. Finally, we note that the transition 



between the generic UV (3.16) and the boundary case (4.2) is somewhat unclear. It appears 



that the solutions first approach the boundary case before switching to the generic behaviour 
in the UV, the scale at which this occurs presumably being associated with the gaugino mass. 
However, more detailed study of the solutions in this region would be necessary to understand 
this completely. 
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Appendix A Equations of motion 



The equations of motion for the full non-SUSY system (section 3.2) can be obtained either 



from the Einstein, Maxwell, dilaton and Bianchi equations of the ten-dimensional system, or 
from a one-dimensional effective Lagrangian L = T — U , with 



2g' (2ti + k' + 2$') + (g'f 
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128 e2$ { e49 (a ') 2 + ( 6 ') 2 ^ " 8e2{g+H) 

+ 2b! (k 1 + 2$') + (ti) 2 + 2$' (k' + }, 

a 4 e 4 9 ( N 2 + e 4k\ _ 4a 3 be 4g N 2 + 2a 2 e 2 9 ( yp^l 



-2{g+h-$) 

+ e 29 N 2 + Ae 2h N 2 c - 8e 2 ^ 9+h+k ^ + 4e 49+2h - e wk + Ae 2h+Ak ) 

- Aabe 2 °N 2 (> + 4e 2h ) + 8b 2 N 2 e 2 ^ + e 4 °N 2 + l6e 4h N 2 

- \Q e 2 ( 2 9+h+k) _ Q4e 2(g+2h+k) + ^(g+k) + lQ e Hh+k)~ _ 



(A.l) 



In addition to the equations of motion resulting from ( A.l ), there is a Hamiltonian constraint 



T + U = resulting from invariance under reparametrisation of the radial coordinate. 



29 



The equations of motion themselves, setting N c = 1 for simplicity, are 

g » = I e -4 9 -2h U, /A 2 _ 4a 2 e 2 g +4fc _ 4a 2 e 2 9 + + 8a6e 2 9 

8 L 

- e 2 s (6') 2 - 46 2 e 2 s - lQe 49+2h g'ti - 16e 49+2h g' & 

- I6e 4 ° +2h (gf + 32e 29+2h+2k - We 2h+4k - 16e 2h ] (A.2) 

h" = - l -e- 2 v- ih [{a') 2 e 49+2h + a 4 e 29+4k + a 4 e 29 - Aa 3 be 29 + Aa 2 b 2 e 2 3 

- 8 a 2 e 2 f +2/l+2fc + Aa 2 e^ +2h - 2a 2 e 29+4k + 2a 2 e 2 ° 
+ 4 a 2 e 2h+4k + Aa 2 e 2h - 4abe 29 - 8abe 2h + e 2h (b'f 

+ Ab 2 e 2h + I6e 29+4h g'ti + 16e 2ff+4 V$' + 16e 29+4h (fc') 2 

_ 8e 2 3 +2^+2fe + e 2 9 +4fc + e 2 9 j (A 3) 

fc" = l - e -^ h {a^ +ik - a 4 e 49 + 4a 3 be 49 - Aa 2 b 2 e 49 + 8a 2 e 29+2h+4k 

_ 8a 2 e 2 9+ 2h _ 8a 2 e S g+ 2h _ 2a 2 e 4 9+ 4fc _ + l 6a6e 2 9 +2ft 

+ Aabe 49 - 8b 2 e 29+2h - 16e 49+4h g'k' - 16e 4g+ih h'k' 

- I6e 49+4h k'<f>' + e 49+4fc - e 49 + I6e 4h+4k - 16e 4ft ) (A.4) 

$" = I e -4 9 -4h r 4 e 4 9 _ 4a 3 be 4g + ^2 ^g + _ 16abe 2g+2h 

8 

+ 2a 2 e 49 - Aabe 49 + 2 (bf e 29+2h + 8b 2 e 29+2h - 16e 49+4h g' & 

- 16e 49+4/ W - lee 49 ^ ($') 2 + e 4g + 16e 4 H (A.5) 

a" = e- A9 - 2h (-\de^ 2h g' - 2a' e 49+2h & + a 3 e 29+4k + a 3 e 29 - 3a 2 be 29 
+ 2ab 2 e 29 - 8ae 29+2h+2k + Aae 49+2h - ae 29+4k + ae 29 
+ 4ae 2?i+4fc + 4ae 2h - be 29 - 46e 2/l ) (A.6) 

b" = -e- 2h (a 3 e 29 - 2a 2 be 29 + ae 29 + Aae 2h + 2e 2 W - Abe 2h ^j (A.7) 
The case discussed in section 4.2, with v% = 0, is far simpler. After setting a = b = 1 and 



g = k the equations of motion for the remaining three functions are 

k" = 2 — 2e~ 4k - 2h'k' - 2{k') 2 - 2k'<f>', 
h" = e 2k ~ 2h - 2h!k' - 2{h') 2 - 2ti&, 

= 2e" 4fe - 2h'& - 2($') 2 - 2fc'* / , (A.8) 

and the constraint is 

e -ik _ e 2k-2h _ 3 + Qh , k , + 4/j / $ / + 6A ./ $ / + ^/j2 + 3 ^/)2 + 2(^')2 = o. (A.9) 

Appendix B Obtaining the UV parameters 



To look for the UV behaviour of the solutions in [To], a matching procedure was proposed 
to provide a fit of the UV parameters. However, we find this process is unreliable when we 
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match at large p. Here we are interested in looking at solutions where the the two scales ph x 
and psusy are varied over a large range. In particular, we need to include cases in which 
one or both have large values. The correct UV behaviour for solutions of this type is only 
manifest at large p, meaning that the matching procedure used in [18J is unsuitable. 

Instead of performing this full matching procedure, it is possible to estimate some param- 
eters from the leading behaviour of appropriate combinations of the background functions. 
For example, we can use the combination 



3 2p/3 a'(p) -> W20 



(B.l) 



to give an approximation of the SUSY breaking parameter in cases where the matching 
procedure fails. Using this method, we find that it appears that W20 — > 00 for vi — > (see 
figure 11). The case W2 = —2 has W^i^) = 0, as would be expected from the invariance of 



W20 under the transformation (4.11). 



Using the same method, the leading coefficients in Mi (see (5.2, 5.3)), can be seen to have 



similar behaviour. However, both these quantities vanish for W2 = 0, with probably non-zero 



derivatives. This reflects the fact that their signs change under the transformation (4.11). 
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Figure 11: Plots of the some of the UV parameters, estimated using the method described 
above for hi = 2.3 with N c = 1 and 0o = (this includes the solutions plotted in figure [8]). 
In (i) we plot W20, corresponding to the gaugino mass which breaks SUSY. We also include 
(ii) = 3W| + 2Q a /c + and (iii) Mf 1 = 8Hn/c+, which contribute to the leading term in 
Mi discussed in section [5] The dotted lines diverging for ii^ — > indicate the position of the 
next point, at w 2 ~ 1.3 x 10" 3 . This has W 20 ~ 10 2 , Mf° w -6 X 10 5 and Mf 1 «3x 10 5 . 



References 

[1] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998), arXiv:hep-th/9711200. 

[2] I. R. Klebanov and M. J. Strassler, JHEP 08, 052 (2000), arXiv:hep-th/0007191. 

[3] J. M. Maldacena and C. Nunez, Phys. Rev. Lett. 86, 588 (2001), arXiv:hep-th/0008001. 

[4] A. H. Chamseddine and M. S. Volkov, Phys. Rev. Lett. 79, 3343 (1997), arXiv:hep- 
th/9707176. 



31 



[5] E. Witten, Adv. Theor. Math. Phys. 2, 505 (1998), arXiv:hep-th/9803131. 
[6] A. Buchel, Nucl. Phys. B600, 219 (2001), arXiv:hep-th/0011146. 

[7] A. Buchel, C. P. Herzog, I. R. Klebanov, L. A. Pando Zayas, and A. A. Tseytlin, JHEP 
04, 033 (2001), arXiv:hep-th/0102105. 

[8] S. S. Gubser, C. P. Herzog, I. R. Klebanov, and A. A. Tseytlin, JHEP 05, 028 (2001), 
arXiv:hep-th/0102172. 

[9] S. S. Gubser, A. A. Tseytlin, and M. S. Volkov, JHEP 09, 017 (2001), arXiv:hep- 
th/0108205. 

[10] N. J. Evans, M. Petrini, and A. Zaffaroni, JHEP 06, 004 (2002), arXiv:hep-th/0203203. 

[11] O. Aharony, E. Schreiber, and J. Sonnenschein, JHEP 04, 011 (2002), arXiv:hep- 
th/0201224. 

[12] R. Apreda (2003), arXiv:hep-th/0301118. 

[13] V. Borokhov and S. S. Gubser, JHEP 05, 034 (2003), arXiv:hep-th/0206098. 

[14] J. Babington, D. E. Crooks, and N. J. Evans, JHEP 02, 024 (2003), arXiv:hep- 
th/0207076. 

[15] J. Babington, D. E. Crooks, and N. J. Evans, Phys. Rev. D67, 066007 (2003), arXiv:hep- 
th/0210068. 

[16] R. Apreda, F. Bigazzi, and A. L. Cotrone, JHEP 12, 042 (2003), arXiv:hep-th/0307055. 
[17] A. Dymarsky and S. Kuperstein (2011), arXiv:1111.1731. 

[18] S. Bennett, E. Caceres, C. Nunez, D. Schofield, and S. Young (2011), arXiv:1111.1727. 

[19] A. Butti, M. Grana, R. Minasian, M. Petrini, and A. Zaffaroni, JHEP 03, 069 (2005), 
arXiv:hep-th/0412187. 

[20] J. Maldacena and D. Martelli, JHEP 01, 104 (2010), arXiv:0906.0591. 

[21] D. Elander, J. Gaillard, C. Nunez, and M. Piai, JHEP 07, 056 (2011), arXiv: 1104.3963. 

[22] R. P. Andrews and N. Dorey, Phys. Lett. B631, 74 (2005), arXiv:hep-th/0505107. 

[23] R. P. Andrews and N. Dorey, Nucl. Phys. B751, 304 (2006), arXiv:hep-th/0601098. 

[24] J. Gaillard, D. Martelli, C. Nunez, and I. Papadimitriou, Nucl. Phys. B843, 1 (2011), 
arXiv: 1004.4638. 

[25] R. Minasian, M. Petrini, and A. Zaffaroni, JHEP 04, 080 (2010), arXiv:0907.5147. 



32 



[26] N. Halmagyi (2010), arXiv:1003.2121. 

[27] E. Caceres, C. Nunez, and L. A. Pando-Zayas, JHEP 03, 054 (2011), arXiv: 1101. 4123. 

[28] R. Casero, C. Nunez, and A. Paredes, Phys. Rev. D73, 086005 (2006), arXiv:hep- 
th/0602027. 

[29] R. Casero, C. Nunez, and A. Paredes, Phys. Rev. D77, 046003 (2008), arXiv:0709.3421. 

[30] C. Hoyos-Badajoz, C. Nunez, and I. Papadimitriou, Phys. Rev. D78, 086005 (2008), 
arXiv:0807.3039. 

[31] C. Nunez, A. Paredes, and A. V. Ramallo, Adv. High Energy Phys. 2010, 196714 (2010), 
arXiv: 1002. 1088. 

[32] F. Bigazzi, A. L. Cotrone, and A. Zaffaroni, Phys. Lett. B519, 269 (2001), arXiv:hep- 
th/0106160. 

[33] J. P. Gauntlett, N. Kim, D. Martelli, and D. Waldram, JHEP 11, 018 (2001), arXiv:hep- 
th/0110034. 

[34] E. G. Gimon, L. A. Pando Zayas, J. Sonnenschein, and M. J. Strassler, JHEP 05, 039 

(2003) , arXiv:hep-th/0212061. 

[35] U. Gursoy and C. Nunez, Nucl. Phys. B725, 45 (2005), arXiv:hep-th/0505100. 

[36] G. Papadopoulos and A. A. Tseytlin, Class. Quant. Grav. 18, 1333 (2001), arXiv:hep- 
th/0012034. 

[37] S. S. Gubser, C. P. Herzog, and I. R. Klebanov, JHEP 09, 036 (2004), arXiv:hep- 
th/0405282. 

[38] S. S. Gubser, C. P. Herzog, and I. R. Klebanov, Comptes Rendus Physique 5, 1031 

(2004) , arXiv:hep-th/0409186. 

[39] A. Dymarsky, I. R. Klebanov, and N. Seiberg, JHEP 01, 155 (2006), arXiv:hep- 
th/0511254. 



33 



